Application of computer simulation free-energy methods to compute the free energy of micellization as a function of micelle composition. 1. Theory

The widespread use of surfactant mixtures and surfactant/solubilizate mixtures in practical applications motivates the development of predictive theoretical approaches to improve fundamental understanding of the behavior of these complex self-assembling systems and to facilitate the design and optimization of new surfactant and surfactant/solubilizate mixtures. This paper is the first of two articles introducing a new computer simulation-free-energy/molecular thermodynamic (CS-FE/MT) model. The two articles explore the application of computer simulation free-energy methods to quantify the thermodynamics associated with mixed surfactant/cosurfactant and surfactant/solubilizate micelle formation in aqueous solution. In this paper (article 1 of the series), a theoretical approach is introduced to use computer simulation free-energy methods to compute the free-energy change associated with changing micelle composition (referred to as DeltaDeltaGi). In this approach, experimental critical micelle concentration (CMC) data, or a molecular thermodynamic model of micelle formation, is first used to evaluate the free energy associated with single (pure) surfactant micelle formation, g(form,single), in which the single surfactant micelle contains only surfactant A molecules. An iterative approach is proposed to combine the estimated value of gform,single with free-energy estimates of DeltaDeltaGi based on computer simulation to determine the optimal free energy of mixed micelle formation, the optimal micelle aggregation number and composition, and the optimal bulk solution composition. After introducing the CS-FE/MT modeling framework, a variety of free-energy methods are briefly reviewed, and the selection of the thermodynamic integration free-energy method is justified and selected to implement the CS-FE/MT model. An alchemical free-energy pathway is proposed to allow evaluation of the free-energy change associated with exchanging a surfactant A molecule with a surfactant/solubilizate B molecule through thermodynamic integration. In article 2 of this series, the implementation of the CS-FE/MT model to make DeltaDeltaGi free-energy predictions for several surfactant/solubilizate systems is discussed, and the predictions of the CS-FE/MT model are compared with the DeltaDeltaGi predictions of a molecular thermodynamic model fitted to relevant experimental data.     

Application of computer simulation free-energy methods to compute the free energy of micellization as a function of micelle composition. 2. Implementation

In this paper, the implementation of the CS-FE/MT model introduced in article 1 is discussed, and computer simulations are performed to evaluate the feasibility of the new theoretical approach. As discussed in article 1, making predictions of surfactant/solubilizate aqueous solution behavior using the CS-FE/MT model requires evaluation of DeltaDeltaG for multiple surfactant-to-solubilizate or surfactant-to-cosurfactant transformations. The central goal of this article is to evaluate the quantitative accuracy of the alchemical computer simulation method used in the CS-FE/MT modeling approach to predict DeltaDeltaG for a single surfactant-to-solubilizate or for a single surfactant-to-cosurfactant transformation. A hybrid single/dual topology approach was used to morph the ionic surfactant sodium dodecyl sulfate (SDS) into the ionic solubilizate ibuprofen (IBU), and a dual topology approach was used to morph the nonionic surfactant octyl glucoside (OG) into the nonionic solubilizate p-aminobenzoate (PAB). In addition, a single topology approach was used to morph the nonionic surfactant n-decyl dimethylphosphine oxide (C10PO) into the nonionic cosurfactant n-decyl methyl sulfoxide (C10SO), the nonionic surfactant octylsulfinyl ethanol (C8SE) into the nonionic cosurfactant decylsulfinyl ethanol (C10SE), and the nonionic surfactant n-decyl methyl sulfoxide (C10SO) into the nonionic cosurfactant n-octyl methyl sulfoxide (C8SO). Each DeltaDeltaG value was computed by using thermodynamic integration to determine the difference in free energy associated with (i) transforming a surfactant molecule of type A into a cosurfactant/solubilizate molecule of type B in a micellar environment (referred to as DeltaG2) and (ii) transforming a surfactant molecule of type A into a cosurfactant/solubilizate molecule of type B in aqueous solution (referred to as DeltaG1). CS-FE/MT model predictions of DeltaDeltaG for each alchemical transformation were made at a number of simulation conditions, including (i) different equilibration times at each value of the coupling parameter lambda, (ii) different data-gathering times at each lambda value, and (iii) simulation at a different number of lambda values. For the three surfactant-to-cosurfactant transformations considered here, the DeltaDeltaG values predicted by the CS-FE/MT model were compared with DeltaDeltaG values predicted by an accurate molecular thermodynamic (MT) model developed by fitting to experimental CMC data. Even after performing lengthy equilibration and data gathering at each lambda value, physically unrealistic values of DeltaDeltaG were predicted by the CS-FE/MT model for the transformations of SDS into IBU and of OG into PAB. However, more physically realistic DeltaDeltaG values were predicted for the transformation of C10PO into C10SO, and reasonable free-energy predictions were obtained for the transformations of C8SE into C10SE and C10SO into C8SO. Each of the surfactant-to-cosurfactant transformations considered here involved less extensive structural changes than the surfactant-to-solubilizate transformations. As computer power increases and as improvements are made to alchemical free-energy methods, it may become possible to apply the CS-FE/MT model to make accurate predictions of the free-energy changes associated with forming multicomponent surfactant and solubilizate micelles in aqueous solution where the chemical structures of the surfactants, cosurfactants, and solubilizates differ significantly.     

Structural, magnetic, high-frequency and high-field EPR investigation of double-stranded heterometallic [{Ni(en)2}2(micro-NCS)4Cd(NCS)2](n).nCH3CN polymer self-assembled from cadmium oxide, nickel thiocyanate and ethylenediamine

A novel heterometallic 1D coordination polymer [{Ni(en)2}2(micro-NCS)4Cd(NCS)2](n) x nCH3CN (en = ethylenediamine) has been prepared using the self-assembly process in a one-pot reaction of cadmium oxide, nickel and ammonium thiocyanates with an acetonitrile solution of ethylenediamine. The complex consists of an uncommon cis-Cd(SCN)4(NCS)2(4-) fragment and a rare combination of cis-Ni(en)2(2+) and trans-Ni(en)2(2+) building blocks linked by micro(1,3)-NCS bridges into a double-stranded zigzag chain structure. Each chain is comprised of [Ni2Cd2(micro-NCS-N,S)4)] macrocycles with chair-like and rectangular-like shapes arrayed alternately. The shortest intrachain CdCd separations are 9.535(1) and 10.868(2) A, while the nearest NiNi distances are 5.418(1) and 6.612(2) A. A network of weak N-HS hydrogen bonds, involving the terminal NCS ligands and NH2-groups of en, links the infinite chains and results in the formation of an extended supramolecular three-dimensional framework. Variable-temperature (1.8-300 K) magnetic susceptibilities show a slight change of the micro(B) value at low temperature, indicative of weak antiferromagnetic interactions (J = 1.55 cm-1) between magnetic canters. High-field, high-frequency (100-400 GHz) EPR spectra were simulated using S = 1 ground state for separate Ni2+ ions with the spin Hamiltonian parameters g = 2.165, D = 0.45 cm-1 and E = 0.03 cm(-1). According to DFT calculations, the D and E parameters are -0.35 cm(-1) and 0.049 cm-1 for the cis arrangement of Ni2+ and 0.58 cm(-1) and 0.012 cm(-1) for trans.     

Structure, formation and catalytic studies of a meso-palladioporphyrin intermediate in a Heck reaction

A meso-palladioporphyrin intermediate was isolated from a Heck reaction between an iodoporphyrin and a non-activated olefin using a Pd(PPh3)2Cl2/Et3N system; its structure was characterized by NMR, MS and X-ray crystallography. Studies on its formation indicate that the Pd(II) catalyst was reduced in situ by Et3N with the assistance of water. The catalytic activity of the intermediate was confirmed by stoichiometric and catalytic reactions using a more reactive olefin, ethyl acrylate.     

Cycloauration of pyridyl sulfonamides

The pyridyl-2-alkylsulfonamides C(5)H(4)N(CH(2))(n)NHSO(2)R (n = 1,2; R = Me, Ph or p-C(6)H(4)Me) and 8-(p-tosylamino)quinoline undergo facile cycloauration reactions with H[AuCl(4)] in water, giving metallacyclic complexes coordinated through the pyridyl (or quinolyl) nitrogen atom and the deprotonated nitrogen of the sulfonamide group. The complexes have been fully characterised by NMR spectroscopy, ESI mass spectrometry and elemental analysis. The X-ray crystal structures of two derivatives reveal the presence of non-planar sulfonamide nitrogen atoms. The complexes show low activity against P388 murine leukaemia cells, possibly as a result of their ease of reduction with mild reducing agents.     

Topical treatment with OGG1 enzyme affects UVB-induced skin carcinogenesis

Nonmelanoma skin cancer resulting from UVB exposure is a large and growing problem in the United States. Production of reactive oxygen species (ROS) during the UVB-induced inflammatory response results in the formation of oxidative DNA adducts such as 8-hydroxy-2-deoxyguanine (8-oxo-dG), which have been shown to contribute to the development of this cancer. The 8-oxoguanine DNA glycosylase (OGG1) enzyme repairs 8-oxo-dG adducts, suggesting that enhancing its activity in the skin might increase 8-oxo-dG repair thus preventing skin cancer development. We therefore used the SKH-1 murine model to examine the effect of topically applied OGG1 on UVB-induced skin cancer development. Mice were exposed three times weekly to UVB followed immediately by topical treatment with a formulation of liposome-encapsulated OGG1 enzyme for 25 weeks. While this treatment did not affect UVB-induced tumor multiplicity, it did reduce tumor size and dramatically reduced tumor progression, as indicated by tumor grade. These results suggest that oxidative DNA damage contributes to the progression of UVB-induced skin tumors and that a topical formulation containing OGG1, perhaps in conjunction with other DNA repair enzymes such as T4 endonuclease V, could be used in populations at high risk for skin cancer development.     

Spectroscopy, binding to liposomes and production of singlet oxygen by porphyrazines with modularly variable water solubility

Three novel classes of porphyrazine-like structures were synthesized to form modular structures in which lipophilicity and water solubility can be tuned. Subtle modification of solubility is an important criterion in selecting a compound for biological photosensitization. The general structure takes the form H2[pz(AnB4-n)], where the core is a porphyrazine (pz) group, A is a pyrrole ring with two sulfide linkages (SR moieties) and B is a pyrrole fused with a 4,7-bis(isopropyloxy)benzo group, with n=4, 3 and 2. These molecules possess their longest wavelength absorption band between 700 and 810 nm, hence laser beams of higher tissue penetration depth could be used to illuminate them in photodynamic therapy (PDT). Armed with absorption bands in the far-red and near-infrared (near-IR), and a capability to tune the solubility, these molecules could make for better sensitizers because of optimized uptake by lipidic membranes and better optical properties. We tested several derivatives of the A4, A3B and A2B2 structures for their singlet oxygen quantum yields in methanol and in liposomes, using 9,10-dimethyl anthracene (DMA) as a singlet oxygen target. Singlet oxygen quantum yields in liposomes ranged from 0.01 to 0.44, with the A2B2 group showing the most promise. In the binding assay to find the equilibrium binding constant, Kb, we detected fluorescence changes due to a change in environment. Peripheral long-chain moieties (the R group in the SR moieties) dominate lipid binding. These moieties range in the hydrophobicity that they induce from C8H17 and benzene, which rendered the molecule totally insoluble in water, to polyethylene glycol (PEG) and carboxylate groups, which imparted water solubility. Each molecule had between 4 and 8 such identical chains. Chains bearing an ether or ester link resulted in measurable equilibrium constants, with a higher Kb for ether substituents. Results for Kb ranged from 0.23 to 26.52 (mg mL(-1))(-1). A delicate balance exists between water solubility and good partitioning to membranes. In general, a higher oxygen-to-carbon ratio in the chains improves binding. Fewer chains and a centrally coordinated zinc ion further improve binding and singlet oxygen production.     

A behavioral model for estimating population exposure to solar ultraviolet radiation

Determining the variability of solar UV exposure of different members of a population by direct measurement demands high compliance over an extended period of time by a large number of people. An alternative approach is to model the variables that affect personal exposure and this is the basis of the method reported here, which uses a random sampling technique to explore variability of exposure at different times of the year by habitués. It is shown that there are large variations in daily personal erythemal exposure, more so for indoor workers living in northern Europe than those resident in Florida, which are due not only to seasonal changes in ambient, but just as importantly to seasonal variation in behavior. Not surprisingly, holiday and summer weekend exposure account for the largest daily UV doses. Northern Europeans who take their summer vacation in Florida can double their exposure during this period compared with holidaying at home and this illustrates just how important sun protection measures should be during recreational exposure in areas of high insolation if the annual UV burden is to be sensibly controlled.